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Photoproduction of charm near threshold 
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Charm and bottom production near threshold is sensitive to the multi-quark, gluonic, and hidden- 
color correlations of hadronic and nuclear wavefunctions in QCD since all of the target's constituents 
must act coherently within the small interaction volume of the heavy quark production subprocess. 
Although such multi-parton subprocess cross sections are suppressed by powers of l/rng, they have 
less phase-space suppression and can dominate the contributions of the leading-twist single-gluon 
subprocesses in the threshold regime. The small rates for open and hidden charm photoproduction 
at threshold call for a dedicated facility. 



PACS: 13.60.Le, 13.60.-r, 12.40. Nn, 12.40. Lg 

The threshold regime of charmonium and open charm 
production can provide a new window into multi-quark, 
gluonic, and hidden-color correlations of hadronic and 
nuclear wavefunctions in QCD. For example, consider 
charm photoproduction 7 p — * J /if) p at the threshold 
energy Elf = 8.20 GeV. [See Fig. @.] The available pro- 
duction energy cannot be wasted at threshold, so all three 
valence quarks of the target nucleon must interact coher- 
ently within the small interaction volume of the heavy 
quark production subprocess. In the case of threshold 
charm photoproduction on a deuteron 7 d — > J/ip d, 
all color configurations of the six valence quarks will be 
involved at the short-distance scale l/m c . Thus the ex- 
changed gluons can couple to a color-octet quark cluster 
and reveal the "hidden-color" part of the nuclear wave 
function, a domain of short-range nuclear physics where 
nucleons lose their identity |l| |j. 

At high energies the dominant contribution to an inclu- 
sive process involving a hard scale Q comes from "leading 
twist" diagrams, characterized by only one parton from 
each colliding particle participating in the large momen- 
tum subprocess. Since the transverse size scale of the 
hard collision is 1/Q, only partons within this distance 
can affect the process. The likelihood that two partons of 
the incident hadrons can be found so close to each other 
is typically proportional to the transverse area 1/Q 2 and 
leads to the suppression of higher- twist, multi-parton 
contributions. However, in contrast to charm produc- 
tion at high energy, charm production near threshold re- 
quires all of the target's constituents to act coherently in 
the heavy quark production process: only compact pro- 
ton Fock states with a radius of order of the Compton 
wavelength of the heavy quark can contribute to charm 
production at threshold. Although the higher-twist sub- 
process cross sections are suppressed by powers of 1/m 2 ., 



they have much less phase-space suppression at thresh- 
old. Thus charm production at threshold is sensitive to 
short-range correlations between the valence quarks of 
the target, and higher-twist multi-gluon exchange reac- 
tions can dominate over the contributions of the leading- 
twist single-gluon subprocesses. 
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FIG. 1. The characteristic scales in elastic J/ip produc- 



tion on protons near threshold, El 



8.20 GeV. The lon- 



gitudinal coherence length of the cc fluctuation of the pho- 
ton is short, l c £S 2E]f/Am 2 c = 0.36 fm. The large mass 
of the charmed quark also imposes a small transverse size 
r± ~ l/m c = 0.13 fm on this fluctuation. The minimum 
momentum transfer is large, t m in ~ —1.7 GeV 2 . All of the 
partons of the target wavefunction have to transfer their en- 
ergy to the charm quarks within their proper creation time 
l/m c , and must be within this transverse distance from the cc 
and from each other, so that charm production near threshold 
occurs at small impact distances 6 ~ l/m c ~ 0.13 fm. 

One can determine the power-law dependence of multi- 
parton heavy quark production subprocesses using an 
operator product analysis of the effective heavy quark 
theory. The heavy quark photoproduction cross section 
can be computed through the optical theorem from the 
corresponding cut diagrams of the forward Compton am- 
plitude. Such diagrams factorize into the convolution of 



two factors: a heavy quark loop diagram connecting the 
photons to the exchanged gluons, times the gauge invari- 
ant matrix element of a product of gluon field strengths 
< p\G™ \p >. Because of the non-Abelian coupling, a sin- 
gle field strength can correspond to one or two exchanged 
gluons. For heavy quark masses, itiq 3> ^qcd ^ ne neav y 
quark loop contracts to an effective local operator, so that 
the field strengths in the matrix element are all evaluated 
at the same local point. The minimal gluon exchange 
contribution (n — 2) gives the leading twist photon- 
gluon fusion contribution. Since < p|G n |p > scales as 
(Ag C£) ) n , each extra gluon field strength connecting 
to the heavy quark loop must give a factor of (l/m|). 
(Higher derivatives in the matrix element are further sup- 
pressed.) Thus one pays a penalty of a factor (A 2 /toq) as 
the number of exchanged gluon fields is increased. How- 
ever, as we shall see, the suppression from the multiple 
gluon exchange contributions are systematically compen- 
sated by fewer powers of energy threshold factors, so that 
at threshold multi-gluon contributions will dominate. A 
similar effective field theory operator analysis has been 
used Q| to estimate the momentum fraction carried by 
intrinsic heavy quarks in the proton |5|||. 

In this paper, we will use reasonable conjectures for 
the short distance behavior of hadronic matter inferred 
from properties of perturbative QCD and effective heavy 
quark field theory to estimate the behavior of the reaction 
cross section. 

The effective proton radius in charm photoproduction 
near threshold can be determined from the following ar- 
gument fj]J|. As indicated in Fig. ||a, most of the pro- 
ton momentum may first be transferred to one (valence) 
quark, followed by a hard subprocess jq — * ccq. If the 
photon energy is E 1 — £E*, where E^ h is the energy 
at kinematic threshold (C > 1), the valence quark must 
carry a fraction x — l/£ of the proton (light-cone) mo- 
mentum. The lifetime of such a Fock state (in the light- 
cone or infinite momentum frame) is r = 1/AE, where 



AE = 1 
2p 
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(1) 



For x = 1/C close to unity such a short lived fluctuation 
can be created (as indicated in Fig. ^a) through momen- 
tum transfers from valence states (where the momentum 
is divided evenly) having commensurate lifetimes r and 
transverse extension 
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(2) 



This effective proton size thus decreases towards thresh- 
old (£ — > 1), reaching r\ ~ 1/wz 2 at threshold (£ — 1 ~ 

As the lifetimes of the contributing Fock states ap- 
proach the time scale of the cc creation process, the time 



ordering of the gluon exchanges implied by Fig. |^a ceases 
to dominate higher-twist contributions such as that of 
Fig. |b @. There are in fact reasons to expect that the 
latter diagrams give a dominant contribution to charmo- 
nium production near threshold. First, there are many 
more such diagrams. Second, they allow the final state 
proton to have a small transverse momentum (the glu- 
ons need p± ~ m c to couple effectively to the cc pair, yet 
the overall transfer can still be small in Fig. ||b). Third, 
with several gluons coupling to the charm quark pair its 
quantum numbers can match those of a given charmo- 
nium state without extra gluon emission. 





(b) 

FIG. 2. Two mechanisms for transferring most of the 
proton momentum to the charm quark pair in 7p — > cap near 
threshold. The leading twist contribution (a) dominates at 
high energies, but becomes comparable to the higher-twist 
contribution (b) close to threshold. 

The above discussion is generic, and does not indicate 
how close to threshold the new effects actually manifest 
themselves. While this question can only be settled by 
experiment, we rely on a simple model to get an estimate 
of the cross section. 

Near-threshold charm production probes the x ~ 1 
configuration in the target, the spectator partons car- 
rying a vanishing fraction x ~ of the target momen- 
tum. This implies that the production rate behaves near 
x — > 1 as (1 — x) 2 " 3 where n s is the number of specta- 
tors |J. Perturbative QCD predicts three different glu- 
onic components of the photoproduction cross-section: 

i) The leading twist (1 — a:) 4 distribution for the process 
717 — > ccq, which leaves two quarks spectators (Fig. ^|a); 

ii) Scattering on two quarks in the proton with a net 

distribution |s , ~/qq — > ccqq, leaving one quark spec- 
tator; iii) Scattering on three quark cluster (Fig. ||b) in 

the proton with a net distribution fe^i , "jqqq ~ * ccqqq, 
leaving no quark spectators. There is some arbitrariness 
in the definition of x close to threshold. We shall use 
x = (2m p M + M 2 )/{s - m 2 ), where s = and M 

is the mass of the cc pair, which has the property x = 1 
at threshold. The relative weight of scattering from mul- 
tiple quarks is given by the probability \/Ft?M. 2 that a 
quark in the proton of radius R ~ 1 fm is found within 
a transverse distance 1 / M (see Ref. JhJ ) . 

The two-gluon exchange contribution produces odd 
C quarkonium jgg — > J/ip, thus permitting exclusive 
7P — > J/tpP production. The photon three-gluon cou- 
pling jggg — ► cc produces a roughly constant term at 



threshold in a/v, where it is expected to dominate (here more specific, x = 0.82 at E™° = 10 GeV and x = 0.69 
v = 1/16tt(s - m 2 p ) 2 is the usual phase space factor). It a t Elf 3 = 12 GeV. So we expect that our model still 

makes sense up to the lowest energy range where exper- 
imental data exist. At higher energies one has to rely on 
the variation of the gluon distribution in the vicinity of 
x ~ to reproduce the steep rise of charm photoproduc- 
tion HQ above E}?° « 100 GeV (x < 0.082). 

As shown in Fig. p|, the threshold dependence of our 



produces the r\ c p, \cP and other C even resonances, but 
also J/-0. 

For elastic charm production (when the proton target 
remains bound) , it is also necessary to take into account 
the recombination of the three valence quarks into the 
proton via its form factor, as well as the coupling of the 
photon to the cc pair. For two gluon exchange the cross 
section of the jp — » J/ ipp takes the form: 



da (1-xf 



2^2 



Fl(t)(s-mi) 



while for three gluon exchange it takes the form: 



(l-x)° 



%(t)(s-m 2 p f 



(3) 



(4) 



where F% g (b) and F 3g (t) are proton form factors that take 
into account the fact that the three target quarks recom- 
bine into the final proton after the emission of two or 
three gluons. While they are analogous to the proton 
elastic form factor Fi(i), they are not known. In the 
numerical applications, we have parameterized them as 
F 2 = exp(1.13i), according to the experimental t de- 
pendency of the cross section [[ll|. The (s — m 2 ) 2 term 
comes from the coupling of the incoming photon to the 
cc pair and the spin-1 nature of gluon exchange (see, 
for instance, Ref. [^2[ ). It compensates the same term 
in the phase space v. The normalization coefficient AT 
is determined assuming that each channel saturates the 
experimental cross section measured at SLAC |l3| and 
Cornell El around E~ = 12 GeV. 



conjectured cross sections U 
the scarce existing data pi 



) and (||) is consistent with 
, |i"3| ]. Indeed, there is also 



evidence [L4| that the energy dependence of the J/ip 
elastic photoproduction cross section at forward angles 
is roughly flat up to E 1 k, 12 GeV, in contrast to the 
steep variation observed at higher energies. More accu- 
rate measurements of the J/ip elastic photoproduction 
cross section up to about 20 GeV are clearly needed. 

The existence of five-quark resonances near threshold 
in the •yp — > pec process jl5| would modify our picture. 
However, the qualitative features of the two- and three- 
gluon-exchange cross sections (which differ by orders of 
magnitude near threshold) should remain valid. 

On few body targets, each exchanged gluon may cou- 
ple to a colored quark cluster and reveal the hidden-color 
part of the nuclear wave function, a domain of short- 
range nuclear physics where nucleons lose their identity. 
The existence of such hidden-color configurations is pre- 
dicted by QCD evolution equations || . It is striking that 
in 7G? — > J/ippn, (Fig. |]), the \B$B% > hidden-color state 
of the deuteron couples so naturally via two gluons to 
the J/ippn final state |lS| ], since the coupling of a single 
gluon to a three-quark cluster turns it from a color octet 
to a singlet. 
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FIG. 3. Variation of the J/ijj photoproduction cross sec- 
tion near threshold. Solid line: two gluon exchange (Eqs. |^). 
Dashed line: three gluon exchange (Eq. ^|). 

Notice that expressions (||) and (f|) are valid in a lim- 
ited energy range near threshold, where x ~ 1. To be 
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FIG. 4. The simplest diagram which reveals a hid- 
den-color state in deuterium Jl8| . 

When the nucleon is embedded in a nuclear medium, 
two mechanisms govern the photo- and electroproduc- 
tion of J/tp mesons. The first, the quasi- free production 
mechanism, contributes the following cross section to the 
jd — > J/ippn reaction, when integrated over the angles 
of the spectator neutron ]19|: 



da 



dtd \ n 



da 
It 



4irn 2 p(\ n |) 



p(\ n \)dn = 1 



(5) 



(0) 



where n | is the momentum of the spectator neutron. 
The nucleon momentum distribution p(\ ft |) in deu- 
terium decreases very quickly pi| with increasing neu- 
tron momentum. Consequently, by selecting high neu- 
tron momenta one can suppress quasi-free production 
and measure inherently nuclear effects. The quasi-free 
contribution in Fig. |^ has been computed with the Paris 
wave function |2Gl of the deuterium. 



7 D — > p n J/t// 
t = - 1 . GeV 2 
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FIG. 5. The variation of the cross-section of the reaction 
yD — > puj/ip against the neutron momentum \n\, at fixed t. 
Solid line: quasi-free contribution. Dashed line: contribution 
of a hidden-color component when its probability is 0.1%. 
Dash-dotted curve: the same for a probability of 1%. 

The second contribution comes from coupling the two 
gluons to separate color octet 3-quark clusters. This con- 
tribution is expected to have a flatter momentum distri- 
bution, since the recoil momentum is shared between two 
nucleons. The corresponding cross section can be roughly 
estimated as : 



da da 
dtd n I dt 



Ann 



*yp—>J/ipp 



Jf(l) 

F?(t) 



(7) 



where the fourth power of the nucleon form factor comes 
from the fact that two nucleons have to recombine, each 
at the momentum transfer t/4 (T^j2^]. We assume that 
the form factor of the transition between a colored cluster 
and the nucleon does not differ too much from the nu- 
cleon form factor |^| and that the recoil momentum is 
equally shared between the two colored clusters whose 
relative wave function is (p cc (^). This component of 
the deuterium wave function has not been measured and 
few predictions are available. As an example and to set 
the order of magnitude, the hidden-color contribution in 
Fig. ^| has been obtained using the Fourier transform of 
the wave function depicted in Fig. 11 of Ref. p3[. Since 



it exhibits a node around 500 MeV/c, a node appears in 
the cross-section around n ~ 1 GeV/c. In a more elab- 
orate calculation the sum over the nucleon internal mo- 
mentum would wash out this node. Anyway, this rough 
estimate shows that the hidden-color component con- 
tribution dominates the cross-section above 0.5 GeV/c. 
The calculation reported in [^3| predicts a probability of 
finding a hidden-color component in the deuterium wave 
function of the order of 0.1%. Fig. || also shows what one 
may expect for a probability around 1%. 

Scattering on colored clusters may dominate sub- 
threshold production, since the high momentum of the 
struck nucleon suppresses the quasi free mechanism. On 
deuterium the threshold for J ftp production is ~ 5.65 
GeV, while on heavy nuclei the threshold is simply the 
J/ij) mass 3.1 GeV. 

Let us close this note with two remarks. At threshold, 
the formation length p5| , p6] | (during which the cc pair 
evolves into a J/ip, after its interaction with a nucleon) 



h 



E 



J/1> 



2m c 



0.22 fm E 7 /GeV (8) 



is around 1 fm, considerably smaller than the size of a 
large nucleus. It is thus possible to determine the scat- 
tering cross section of a full sized charmonium on a nu- 
cleon using nuclear targets, in contrast to the situation 
at higher energies where the nuclear interaction of a com- 
pact cc pair is measured. The study of the A dependence 
of the J /ip photoproduction cross section at SLAC at 20 
GeV |24j gave aj/^, N = 3.5 ± 0.8 ± 0.5 mb. However, a 
large calculated background was subtracted and the lack 
of information on the J /ip kinematics prevented a sepa- 
rate measurement of coherent and incoherent photopro- 
duction. A new measurement of J j 1 ^ photoproduction on 
several nuclei around E 1 ~ 10 GeV, with good particle 
identification and a determination of the J/ip momentum 
is clearly called for. 

Although the cc pair is created with rather high mo- 
mentum at threshold, it may be possible to observe re- 
actions where the pair is captured by the target nucleus, 
forming "nuclear-bound quarkonium" JlB") . This process 
should be enhanced in subthreshold reactions. There is 
no Pauli blocking for charm quarks in nuclei, and it has 
been estimated that there is a large attractive Van der 
Waals potential binding the pair to the nucleus |27]] . The 
discovery of such qualitatively new states of matter would 
be very important. 

In this paper we have shown that charm production 
near threshold has strong sensitivity to the multi-quark, 
gluonic, and hidden-color correlations of hadronic and 
nuclear wavefunctions in QCD. Although multi-parton 
subprocess cross sections are suppressed by powers of 
1/m^, they have correspondingly less phase-space sup- 
pression and thus can dominate the contributions of the 
leading-twist single-gluon subprocesses in the threshold 



regime. Such processes will therefore add to our under- 
standing of the short-range structure of nucleons and 
nuclei. A new dedicated facility, such as CEBAF12 or 
ELFE, with high intensity and duty factor, will make 
possible the experiments needed to explore this QCD 
frontier. 
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